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o 

We present a novel isospectral deformation chain in supersymmetric quantum 
mechanics. The chain is applied to the Coulomb potential. 

<: 

I. INTRODUCTION 

n 

d- In 111 we introduced within the framework of supersymmetric quantum mechanics an 

isospectral deformation on the form 

W(x) -> W (x) = F (x)W(x) , (1) 

where W(x) is some known superpotential and F (x) some function to be determined by 
00 ■ the isospectrality condition 

o; 

V + (x) = W%{x) + W' Q {x) = W 2 {x) + W\x) EE V + (x) . (2) 

It was shown that Eq.(l) includes the only previously explored deformation of this kind, 
which has the form ([2j, e.g.) 

W{x) ^W (x) = W(x) + f(x). (3) 

f(x) is some function which is determined by Eq.(2). In this paper we expand the deforma- 
tion Eq.(l) in various directions and study the implications drawn from the isospectrality 
condition. We show in particular that the deformation Eq.(l) is the root of an infinitely 
long and recursively generated chain of deformations. 

In the next section we briefly review some findings in [l|. We then define and derive 
an explicit recursive scheme for generating novel isospectral deformations. We then discuss 
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possible generalizations of this scheme. We next apply the recursive scheme to the Coulomb 
potential. We conclude and discuss some of our findings in the last section. 



II. BASE DEFORMATIONS 



The deformation Eq.(l) implies the following differential equation for F ( 

4-F (x) + (-^-\nW(x))F (x) + W(x)F*(x) = W(x) + ^-\nW(x) . (4) 

ax ax ax 

This is the generalized Riccati equation . If one particular solution F 00 (x) of Eq.(4) is 
known another solution is given by J] 

F (x) = Foo(x) + -i- , (5) 

-Ao{X) 

where X (x) solves the equation 

-^X (x) - In W(x) + 2F 00 (x)W(x))X (x) = W(x) . (6) 
q.(6) can be solved by elementary means. The resulting superpotential Wq(x) is given by 



W Q (x) = (Foo(x) + ■^ n )W(x) = W 00 (x) + y^W(x) = 
A.o{x) ^-o{x) 

e -2f x Fao(t)W(t)dt 

Cqi is an integration constant, which we will assume to be real. We have explicitly introduced 
upper integration limits in Eq.(7) in order to avoid sign ambiguities. This explains the 
difference in the sign in the denominator in Eq.(7) compared with (2.5) in [l| where the 
reverse order of integration in one of the integrals was implicitly assumed. We do not 
specify the lower integration limits in Eq.(7). These are not important, of course, since the 
values of the integrals there can essentially be absorbed into Cq\. When Fq Q (x) = 1 we 
identically rederive Eq.(3) and the corresponding expression discussed in ^|. We can by 
simple inspection see that the particular solution F (x) = 1, the identity deformation, solves 



1 We will often rewrite fractions on the form W'(x)/W(x) as the logarithmic derivative of W(x) as a 
formal tool. Caution must of course be exercised when using the corresponding expressions in actual 
computations. 
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Eq.(4). This deformation corresponds to the limit C i — > oo in Eq.(7) with F 00 (x) = 1. 
In the limit C i — > oo we generally get W (x) = Wo (x). This deformation will play a 
pivotal role in this work; it will represent the base of a recursive scheme for generating novel 
isospectral deformations. We will therefore refer to a particular Woo(x) as a base deformation 
in the following. 

In order to expand the space of concrete isospectral deformations further we transform 
Eq.(4) into an ordinary second order differential equation by the substitution 

This substitution gives rise to the following linear homogeneous equation 

--^U (x) + V+(x)U (x) = 0. (9) 

The special solution F (x) = 1 is generated by the solution 

[/ (x)~e^^. (10) 

The particular solutions for Fq(x) stemming from these equations can be fed into Eq.(7) 
(as F o(x)) and thus expand the space of available concrete deformations. The physical 
potential V-(x) generated by W (x) can in general thus be written [l| 2 

Vt(x) = W*(x) - W&x) = W 2 m {x) - WUx) + 



AW 00 (x)e- 2 I xWo ^ dt 
+ C m + j x e- 2 rw 00 {t)dt du 

with 



-2f x W 00 (t)dt 



n 2 



C m + j x e~ 2 rw 00 (t)dt du 



(11) 



W 00 (x) = ^-\nU (x) . (12) 



III. RECURSIVE LINEAR DEFORMATIONS 

Although the Riccati equation can be transformed into an ordinary second order differ- 
ential equation the non-linearity of the equation allows for a solution space which is larger 
than the one associated with linear differential equations of second order, as became evident 
in the previous section. It is therefore natural to ask whether the non-linearity of the Riccati 
equation implies even more isospectral deformations than the ones we already have deduced 
l|. We will explore this question in this and the next section. 



Note that the corresponding expression in [l| ((2.14)) is misprinted. 
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A. The sum 

Let us entertain the following idea. Assume that we have derived a particular base 
deformation Wqo(x) from an explicitly given superpotential W(x). Then assume that we 
add another term F\{x)W{x) to that deformation such that we in principle get a novel 
deformation on the form W(x) = F w (x)W (x) + Wqq(x). After determining F w (x) from the 
isospectrality condition Eq. (2) add yet another term of this kind to the deformation. Let us 
assume that this process can be repeated indefinitely. Will terms added in this manner give 
rise to novel deformations? We will in the following show that they do. This represents a 
recursive deformation scheme. 

Following the basic idea, after m iterations we thus have the general recursive linear (in 
W(x)) deformation 

m 

W m o(x) = (J2 \F i0 (x))W(x) = X m F m0 (x)W(x) + W {m _ 1)0 (x) , A = 1 . (13) 

i=0 

The Aj's are assumed to be independent real constants. Starting with a known superpotential 
m consecutive applications of the isospectrality condition yields the following set of equations 

F^(x) + [\nW(x)]'F 00 (x) + W{x)F^{x) = W(x) + (ln^(i))' , 
F[ (x) + [(In W(x))' + 2F 00 (x)W(x)}F 10 (x) + X 1 W(x)Ff (x) = , 
F> (x) + [(In W(x)Y + 2(F 00 (x) + A 1 F 10 (x))W(x)]F 20 (x) + \ 2 W(x)F? (x) = , 

(14) 



[ FU*) + [(In W(x))> + 2W {m _ 1)0 (x))]F m0 (x) + X m W(x)F^ (x) = . 
The first equation in Eq.(14) coincides of course per definition with Eq.(4). Note that 
Fj (x) = 1 only solves the first equation in Eq.(14). Let us consider an arbitrary iteration 
level n (7^ 0) and make the following substitution in Eq.(14) 

F n0 (x) = ^(\nU n (x))' . (15) 

The equation for F n (x) can then be written 

K(x) + (X n - 1)^^ + 2W {n _ 1)0 (x)U' n (x) = . (16) 

This equation corresponds to Eq.(9) in the case when n = 0. It reduces in general to an 
ordinary linear differential equation only when \ n = 1 , Wn 7^ 0. We will focus on this special 
case in this paper. 
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The general solution of Eq.(16) for arbitrary n / 0, and with \ n set to unity, can be 
found by elementary means, and we deduce that 

C 2 e~ 2 f x w ( „-i) (t)dt 



C n \ and C n2 are integration constants, which we assume to be real. We can reduce the 
number of integration constants to one at each iteration level, but we will stick to the habit of 
explicitly writing down the actual number of constants in order to make it easier to compare 
the various formulas we deduce, which stem from both second and first order differential 
equations. We also note that the structure of F n0 (x) implies that previous deformations are 
not regenerated in general. Of course, this does not exclude this possibility to arise, as we 
will see in Section 5. Hence, m in Eq.(13) has in principle no natural upper bound. From 
Eq.(13) and Eq.(17) we get the following expression for the superpotential at iteration level 
m 



F n0 (x)W(x) 



Cm + C n2 J x e~ 2 f u ^(»-i)o(0*d u 




(17) 





W Q0 (x) + —lnP m (x). 



(18) 



From Eq.(18) we deduce that 



(19) 



such that 



Pi Or) 
P 2 (x) 



Cn + C 12 J x e- 2 f u ^ dt du, 

Pi(x)(C 21 + C 22 J x P{ 2 {u)e~ 2 r ^ dt du) , 



< 



(20) 



Pm{x) 



P m -i(x)(C m i + C m2 J 



P m \{u)e- 2 rw^du). 



Hence, 



Pn(x) = n(C,-i + C j2 / P--\(u)e- 2 ^ ^ dt du) . 




(21) 
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ttt xF =F o v T J 7 +F 10 W T i, T i, D v +F20W +F m0 W T j, 

VF > Woo > Wio = Woo + (In Pi) > ■ ■ ■ > W m0 



W > Woo = W > Wio = W + (In Pi) s- ••• > W m o 

FIG. 1: The upper line depicts the solvable deformation chain Eq.(13) to iteration level m. There 
is no upper bound on m. The Fjo(x) functions are given in Eq.(17). The Pj(x) functions are given 
in Eq.(18) and Eq.(20). They are functions of a base deformation Woo(x). A base deformation 
Woo(x) is generated by the zero-energy Schrodinger equation interacting with the partner potential 
V+(x), Eq.(9). The second line depicts the important special case when Fo(x) = 1. This particular 
solution can be derived as a special case of Eq.(7) with Xq 1 = 0, which can be achieved by 
Cqi — > 00, and Foo{x) determined by Eq.(8) and the solution Eq.(10). Note that Wio{x) then 
coincides with Eq.(7) when Foo(x) = 1 and Cqi is finite in Eq.(7); Eq.(7) is thus in this particular 
case regenerated by the scheme at the next recursion level, i.e. 

This last form of the P n (x) functions neatly exhibits how the base deformation W o(x) 
generates the higher order deformations. Some of the details we have deduced so far in this 
subsection are presented in Figure 1. 

Make the following substitution at each iteration level in Eq.(13) 

F n o{x) -> F n o{x) + — ^— . (22) 
Ji n yx) 

This implies (with the A m 's reinstated in Eq.(13)) a generalized form W m (x) of the super- 
potential W m o(x) 

m m > 

W m0 (x) = KF i0 W(x) W m {x) = W m0 (x) + YT~) W ^ X) ■ (23) 

From Eq.(14) we find that X n {x) satisfies the equation 

^X n {x) - {j- \nW{x) + 2W n0 (x))X n (x) = \ m W{x) . (24) 

This equation is a generalization of Eq.(6). The n'th deformation term Eq.(17) thus changes 
into 

F n0 (x)W(x) -> F n0 (x)W(x) + ^ln(C n3 + K f e~ 2 ^ W ^ dt du) . (25) 
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C n3 are integrations constants, which we assume to be real. Eq.(25) implies that the more 
general expression for the superpotential in Eq.(23) can be written as 

W m {x) = W m0 (x) + A-\nQ m (x) , (26) 

where 

Qrn{x) = f[(C i3 + Xi f e- 2 rW M {t) d t du) (27) 

i=o 

m = in Eq.(26) (A = 1) reproduces Eq.(7). In the special case when \ m = l,Vm in 
Eq.(13) we get 

W m (x) = Woo(x) + -^\nP m (x) + -^-\nQ m (x) . (28) 

When we compare the expressions for (P m (x))' and (Q m (x))' we find that they differ by just 
the last term in (Q m (x))' '. 



B. The product 

What happens if we in Eq.(13) assume a product structure instead of a sum structure ? 
Let us assume that we have determined a base deformation. Let this be the seed superpo- 
tential for the deformation 

Woo(x) -)■ W 10 (x) = F 10 (x)W 00 (x) = F 10 (x)F 00 (x)W(x) , (29) 

where F w (x) is some function to be determined by the isospectrality condition. This product 
scheme can of course in principle be repeated an arbitrary number m times 

m 

W m0 (x) = (]jF i0 (x))W(x) = F m0 (x)W {m - 1)0 . (30) 

i=0 

This structure gives rise to the following set of equations 

' F^ (x) + (\nW(x))'F 00 (x) + W(x)F* (x) = W(x) + (\nW(x))' , 

F{ (x) + (\nW o(x))'F 10 (x) + W 00 (x)F? (x) = -^(W + (\nW(x))') , 
< : 

FUx) + (lnW {m _ 1)0 (x))'F m0 (x) + W {m _ 1)0 (x)FUx) = w m ^ (x) (W(x) + (\nW(x))') . 

(31) 
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Clearly, each iteration level depends on all the previous ones, and at each level we are 
dealing with a non-homogenous non-linear differential equation. Interestingly, by making 
the following substitution at an arbitrary iteration level n^O 

Fno(x)= , 1 - (hxU n (x))', (32) 

tV(n-l)0W 

where U n (x) is some function, the equations Eq.(31) all reduce to Eq.(9). Hence, attempting 
to generate novel deformations recursively via a product structure, of the kind above, fails. 
This conclusion was also reached in [l| , but at the level of the second order linear differential 
equation Eq.(9). 



IV. RECURSIVE NON-LINEAR DEFORMATIONS 

We have so far only considered linear deformations. In this section we will briefly consider 
two non-linear deformation schemes. Let us first consider a polynomial kind of deformation. 
That is, given a superpotential W(i-i)o(x) which we will assume is derived, in some way or 
another, from some seed superpotential W(x). Consider then the polynomial deformation 

W i0 (x) = F l0 (x)W\x) + W (i -i )0 (x) ; e {1,2,3,...}. (33) 

The isospectrality condition implies 

F; (x) + [k(\n W{x))' + 2%i )0 (z)]F«)(z) + W k {x)F%{x) = . (34) 

This is a Riccati type equation of the kind we have met earlier in this paper. Apparently, 
different fc-values give rise to very different equations to solve. However, and rather intrigu- 
ingly, all the possible fc-values implies the same deformation. This is seen by making the 
following substitution 

W h (x) Uk\x) 

where Uk(x) is some function. This expression inserted into Eq.(34) gives 

U%(x) + 2Wi^{x)U' k {x) = . (36) 

Hence, Wio(x) is independent of k and we are essentially left with a linear deformation. 
Clearly, the range of values of k can be expanded to the real numbers. 
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Another canonical generalization of our work is to consider deformations on the form 

W(x) = H (x)T(W) , (37) 

where T is any functional of the seed superpotential W[x). The isospectrality condition 
then implies 

H' (x) + QnF(W))'H (x) + J r (W)H 2 (x) = F(Wy l (W 2 (x) + W\x)) . (38) 

Note that H (x) = 1 does not solve this equation unless T{W) = W , since Eq.(38) with 
H (x) = 1 implies J-'{x) +J r2 (x) = V+(x). Since V+(x) is uniquely given in Eq.(2) any other 
choice of functional will fail to satisfy the isospectrality condition. Hence the conclusion. 
The particular solution H (x) = 1 is not forced upon us. We can in principle do without it. 
It is easily verified that Eq.(38) can be cast into the form Eq.(9) by the substitution 

H (x) = y^QnU(x))' . (39) 

We can also look for an expanded solution by writing 

H (x) = H 00 (x) + -i-, (40) 
Z (x) 

where H 00 (x) is a particular solution of Eq.(38). We then get the equation 

^-Z (x) - {j- \nT{W) + 2H 00 (x)T(W))Z (x) = JF{W) , (41) 
which is a generalized form of Eq.(6). The reciprocal solution has the general form 

I e -2 f x H 00 (t)T(W)dt 

Z (x) = F(W)(C + J x e-tfHwVrWWdu) ' (42) 

C is an integration constant, which we assume to be real. Utilizing that H Q0 (x) = 
J r ^ 1 (x)(\nU(x))' the resulting deformation coincides with Eq.(7). We thus therefore con- 
clude that non-linear deformations on the form Eq.(37) does not generate additional defor- 
mations to the ones already generated by Eq.(l). 



V. DEFORMING THE COULOMB POTENTIAL 



As a relatively simple application of the linear deformation scheme let us briefly con- 
sider deformations of the Coulomb potential. This potential has, within the framework 
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of supersymmetric quantum mechanics, been treated in several previous works 2|. The 
superpotential and the partner potential for the Coulomb potential are given by {2] 

^ = 2^-^ < 43 > 

q and I in these expressions are the electric charge and the angular momentum quantum 
numbers, respectively. These potentials result in the following general solution for Uq(x) in 
Eq.(9) (J 

U (x) = dM^jC^-) + C 2 W l+1 , +l (^) . (45) 

The M(x)- and W / (x)-functions are the Whittaker functions. The solution Eq.(10) is given 
by fl 

U (x) ~ e*)" (,+1)In(2l) . (46) 

We will for simplicity assume this solution in the following. We will let Cqi — > 00 in Eq.(7) 
such that we deal with the identity deformation Wq(x) = W 00 (x) = W(x). We will also 
ignore the Qj(x) contributions in the following. Define A = q 2 /(2(l + 1)) and B = I + 1. It 
then follows that 

p x ( x ) =C n + C 12 J X t 2B e~ 2At dt , (47) 

such that 

ft ff 2B -2Ax 

(48) 

Let us consider the s-state with I = in order to get a better grasp on the content buried 
in Eq.(48). We also set q = 1. The expression for W / i (x) then reduces to 

T?r . , 1 1 C 12 x 2 e~ x 

w ^ = 2 " s + ^hwP (49) 

after redefining Cn such that the lower integration limit of the integral in Eq.(48) does 
not appear explicitly in the expression for the potential. We will automatically do such 
redefinitions in the following when it is appropriate. The corresponding physical potential 
V-i(x) can either be derived from the definition VLi(x) = W 2 (x) — W{ (x) or from Eq.(ll) 
with Woo(x) = W(x) and Cqi finite. This is a consequence of a regeneration of Eq.(7) by 
the recursion scheme which was noted in Figure 1. From the definition it follows that 

V ( T ) = ---+ ° 12X{2X ~ A)e ~ X + 2C 2 2 x*e- 2 * 
~ 1{ ' 4 x + C n -C 12 (x 2 + 2x + 2)e-* (C n - C 12 (x 2 + 2x + 2)e-) 2 " 1 ) 



11 



In the special case when we set C u = the last term in Eq.(49) becomes independent of the 
exponentials (and C\ 2 ) and thus reduces to a pure rational function. The physical potential 

V-i(x) generated by W\o(x) is then given by 

4x(x + 2) . . 4x(x + 2) 



4 x + {x 2 + 2x + 2) 2 



= VJx) + 



(51) 



(x 2 + 2x + 2) 2 ' 

Let us go to the second iteration level starting from the expression for W w (x) in Eq.(49) 

with Cn = 0, for convenience. It then follows that 

C 22 x 2 e x 



W 20 (x) = W 1 (x) + 



(52) 



C 21 (x 2 + 2x + 2) 2 + C 22 (x 2 + 2x + 2)e* 
Note that when C 2 \ = we get W 2 q(x) = W(x). Hence, the deformation scheme allows 
in general for the possibility that additional iterations in particular cases may regenerate 
previous potentials in a nontrivial fashion. The expression for the corresponding physical 
potential is given by 

C 22 x 2 e x 



V- 2 (x) = V^(x) + 



C 21 (x 2 + 2x + 2) 2 + C 22 (x 2 + 2x + 2)e 



x 



X 



M- + 

y x x 2 + 2x + 2' 



1™2 
2 X 



+ 



2(C 21 (2x + 2) + C 22 e 3 



C 21 (x 2 + 2x + 2) + C 22 e x \ ' ^ 
The superpotential stemming from the third iteration with W 20 (x) in Eq.(52) as the 
starting point is given by 

C 32 x 2 e x 



W 30 (x) = W 20 (x) + 



C 21 (x 2 + 2x + 2) 2 + C 22 (x 2 + 2x + 2)e* 



C 31 + ^ln 



C 21 (x 2 + 2x + 2) + C 22 e* 



(54) 



x 2 + 2x + 2 

This superpotential introduces the possibility for a logarithmic singularity away from the 
origin when C 22 /C 21 < 0. We note that setting C 31 = does not regenerate a previous 
potential as was possible at the previous iteration level when we correspondingly put C 2 \ = 0. 
From Eq.(54) we can deduce the physical potential ^-3(2;) at the third iteration level. We 
do not reproduce it here due to its complexity. Due to the complicated integrals appearing 
we are not able to provide the analytical expression for W4o(x). We leave detailed studies 
of the Coulomb potential for the future. 



VI. CONCLUSION 



In a previous paper we showed that isospectral deformations on the form Eq.(3) are 
contained in the space of deformations generated by isospectral deformations on the form 
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Eq.(l). In this paper we have shown that Eq.(l) can be considered as the initial, or base, 
deformation of a novel infinite recursive isospectral deformation chain. This thus answers 
to some extend the question by which we ended our previous paper Qj; how does the most 
general isospectral deformation of the kind considered there (Eq.(l) in this paper) look 
like. The results in this paper do obviously only give a partial answer. We deduced in 
particular that a class of recursive deformations exists which is generated by the solutions 
of the non-linear differential equations in Eq.(16). 

We briefly discussed various ways to construct alternative recursive deformation struc- 
tures. We considered a linear product structure, polynomial deformations and completely 
generalized base deformations. They all either failed to provide a recursive structure or they 
turned out to be identical to the linear deformation scheme. 

We applied the linear recursive scheme to the Coulomb potential. We derived novel su- 
perpotentials which all satisfy the isospectrality condition. This application did also demon- 
strate how easily novel isospectral deformations can be generated in this approach. It did 
also demonstrate an increased relative complexity of the generated potentials with the num- 
ber of iterations, as one also naively would expect from the expression Eq.(21). 
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